Probing Non-Abelian Statistics with Quasiparticle Interferometry 
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We examine interferometric experiments in systems that exhibit non-Abelian braiding statistics, expressing 
outcomes in terms of the modular S-matrix. In particular, this result applies to FQH interferometry, and we give 
a detailed treatment of the Read-Rezayi states, providing explicit predictions for the recently observed v = 12/5 
plateau. 
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Quantum systems in two spatial dimensions allow for ex- 
otic exchange statistics, characterized by a unitary represen- 
tation of the braid group (jj, |2J. This representation may be 
non-Abelian, acting on a multi-dimensional internal Hilbert 
space [3, 4]. So far, experimental evidence for the existence 
of such anyonic statistics has only recently been found in the 
(Abelian) Laughlin states of Fractional Quantum Hall (FQH) 
systems Qy. However, the prospect of non-Abelian statis- 
tics is far more exciting, especially in light of its potential ap- 
plication in topologically fault- tolerant quantum computation 
00]. There are currently several observed FQH states, at 
filling fractions v = 5/2, 7/2, 12/5 §1 GU EH (and possi- 
bly v — 3/8, 19/8), that are ex pec ted to possess non-Abelian 
statistics. Numerical studies ill 2l Il3l ll4Tl suggest that the 
v = 5/2, 12/5 states should be described respectively by the 
Moore-Read (MR) state and the k = 3, M = 1 Read- 
Rezayi (RR) state 13. Clearly, as experimental capabili- 
ties progress, it becomes increasingly important to understand 
how to probe and correctly identify the braiding statistics of 
quasiparticles. In this Letter, we explain how, for any system 
described by a topological quantum field theory (TQFT) in 
the infrared limit (e.g. FQH systems), knowledge of modular 
S-matrices may be used to extract this information from inter- 
ferometry experiments, or, inversely, to predict the outcomes 
of such experiments. As an example relevant to current exper- 
imental interests, we obtain explicit results for the RR states. 

The topological properties of 2D quantum systems with 
an energy gap can be described using TQFTs (or "modular 
tensor c ateg ories" in mathematicians' terminology, see e.g. 
J3 G3 EJd), often abstracted from conformal field the- 
ories (CFTs, see lEflll and references therein). Such an anyon 
model is defined by (i) a finite set C of particle types or "any- 
onic charges," (ii) fusion rules specifying how these particle 
types may combine or split, and (iii) braiding rules dictat- 
ing the behavior under exchange of two particles (all subject 
to certain consistency conditions). The "vacuum" charge is 
given the label 1. The anti -particle or "charge conjugate" of 
a particle type a is denoted a, and is the unique charge that 
can fuse with a to give 1. Fusion of particle types gener- 
alizes the addition of charges or angular momenta and the 
(commutative and associative) fusion rules are specified as 
ax b = J2cec ^ab c > wnere the integer N£ b is the dimension 
of the Hilbert space of particles of type a and b restricted to 



have total anyonic charge c. Fusion and braiding can be repre- 
sented diagrammatically on oriented, labeled particle world- 
lines, and are unaffected by smooth deformations in which 
the lines do not intersect. Charge conjugation is represented 
by reversal of wordline orientation. We will refer to only one 
braiding relation, known as the modular S-matrix, defined by 
the following diagram: 



Sab = 



*.GD. 



(i) 



Here D — y/J2 a — V^l 1 * s tne total quantum dimen- 
sion, where d a , the quantum dimension of particle type a, is 
the value of a single loop of that type, 



d a = 




DS 



la • 



Some useful properties of the S-matrix are 



Sab — Sba — {S 1 ) a ( ) — S a b 



(2) 



(3) 



The importance of the S-matrix becomes clear when one envi- 
sions interferometry experiments for these systems in which a 
particle has two possible paths that it may take around another 
particle, the two paths combining to form a closed loop. This 
is typical of Mach-Zender, two-s lit, FQH two-point-contact, 
etc. experiments J^H El El |M 111 EE 113 EE El E3l • In 
such experiments, an interference term arises that can be writ- 
ten as 



U^U 2 |* a6 ) = e* 



(4) 

where l^afc) is the initial state of particles a and b, and Ui,Uz 
are the unitary evolution operators for the particle a traveling 
around the particle b via the two respective paths. It has been 
rewritten in terms of the monodromy operator M that contains 
only the contribution from adiabatically transporting particle 
a around particle b (i.e. braiding), and a phase e laah that 
absorbs all other contributions (i.e. it contains the free par- 
ticle dynamics and the Aharonov-Bohm phase from a back- 
ground magnetic flux). For simplicity, we let the two particles 
have definite anyonic charge, however it is a straightforward 
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generalization to allow superpositions of particle type. If the 
theory only has Abelian statistics, then |Af a &| = 1, but with 
non-Abelian statistics, \M a b\ can be less than 1 and must be 
calculated by TQFT methods. The braiding term is diagram- 
matically represented by winding the worldline of particle a 
around that of particle b, taking the standard closure (where 
the worldline of each particle is closed back on itself in a man- 
ner that introduces no additional braiding) and dividing by the 
quantum dimension of the two particle types fi^ll . Thus, the 
resulting monodromy matrix element can be written entirely 
in terms of the S-matrix: 

m «»=tt C(T) & 

This result is particularly nice because the S-matrix is typi- 
cally more readily computable than the complete set of braid- 
ing/fusion rules, and, in fact, has already been computed for 
most physically relevant theories. In particular, this applies to 
the class of theories described by CFTs generated as products 
and cosets of Wess-Zumino-Witten theories, which includes 
all proposed non-Abelian FQH states (see jHH IH 0. H). 
In the case of the product of two CFTs or TQFTs, the par- 
ticle types are denoted by pairs of particle labels, one from 
each theory, and the S-matrix of the product theory is the ten- 
sor product of the S-matrices of the parent theories. In the 
Gk/Hi coset theory, the new labels are also formed as pairs 
of labels from the parent Gk and Hi theories, but now there are 
branching rules which restrict the allowed pairings. Also, dif- 
ferent pairs of labels may sometimes turn out to represent the 
same particle type, a phenomenon described by "field iden- 
tifications". This means that there will be only one row and 
column in the S-matrix for any set of identified labels. De- 
spite these complications, the coset's S-matrix elements are 
described by the simple formula l35ll 

S { S! P fL)=<G,H,k,l)S ! ^S^ (6) 

where a, b and p, q are respectively labels of the Gk and Hi 
theories, and c(G, H, k, I) is an overall normalization con- 
stant that enforces unitarity of the S-matrix (but is irrelevant 
in Eq. Q). 

We now focus specifically on FQH systems, because they 
are the only physical systems found to exhibit braid statistics 
thus far, and represent the most likely candidates for finding 
non-Abelian statistics. We consider the interferometry experi- 
ment originally proposed in lEill for measuring braiding statis- 
tics in the Abelian FQH states, which was later adopted for the 
non-Abelian case in 1 26] and addressed again in the context of 
the 1/ ^ 5/2 state in Il29ll30ll . A somewhat similar experiment 
has recently been implemented to probe the Abelian v = 1/3 
state QEJ]. The experimental setup is a two point-contact in- 
terferometer composed of a quantum Hall bar with two front 
gates on either side of an antidot (see Fig. [0- Biasing the 
front gates, one may create constrictions in the Hall bar, ad- 
justing the tunneling amplitudes t\ and t<2,. Tunneling between 



the opposite edge currents leads to a deviation of a xy from its 
quantized value, or equivalently, to the appearance of a xx . By 
measuring a xx one effectively measures the interference be- 
tween the two tunneling paths around the antidot. The tunnel- 
ing amplitudes t\ and £ 2 must be kept small, to ensure that the 
tunneling current is completely due to quasiholes rather than 
higher charge composites |43], and to allow us to restrict our 
attention to the lowest order winding term. 

In order to be able to influence the resulting interference 
pattern, we envision several experimentally variable parame- 
ters: (i) the central gate voltage allowing one to control the 
number n of quasiholes on the antidot, (ii) the perpendicular 
magnetic field B, (iii) the back gate voltage controlling the 
uniform electron density, and (iv) the side gate that can be 
used to modify the shape of the edge inside the interferomet- 
ric loop I29ll . The intention is to be able to separately affect 
the Abelian Aharonov-Bohm phase and the number of quasi- 
holes on the antidot; from this point of view having all these 
controls is redundant, but may be found beneficial for experi- 
mental success. 




FIG. 1 : A two point-contact interferometer for measuring braiding 
statistics. The hatched region contains an incompressible FQH liq- 
uid. The front gates (F) are used to bring the opposite edge cur- 
rents (indicated by arrows) close to each other to form two tunneling 
junctions. Applying voltage to the central gate creates an antidot in 
the middle and controls the number n of quasiholes contained there. 
An additional side gate (S) can be used to change the shape and the 
length of one of the paths in the interferometer. 

The longitudinal conductivity is proportional to the proba- 
bility that current entering the bottom edge will leave through 
the top edge, which to lowest order is: 

a xx <x \h\ 2 + \t 2 \ 2 + 2Re{t* 1 t 2 {^ ab \U^ 1 U 2 \^ab}} 
- |<i| 2 + |< 2 | 2 + 2|< 1 i 2 ||A/ ofc |cos(/3 + M . (7) 

In this equation (3 = a a b + arg(i 2 /ii) can be varied by 
changing B (keeping the quasihole number fixed), the rela- 
tive tunneling phase, and/or the edge shape around the central 
region. We have written M a f, = \M a b\ e l6ab , and will be in- 
terested in the elements where the particle a which tunnels 
carries the anyonic charge of a quasihole, while the "particle" 
b is a composite of n quasiholes on the antidot, carrying an 
anyonic charge allowed by the fusion rules |44|. 

We will now apply this formalism to the RR states and, 
in particular, look more closely at the k = 3, M = 1 case, 
which, by a particle-hole transformation (which generally in- 
verts the statistics and has the effect of conjugating the S- 
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matrix), is the expected description of v = 12/5. The anyon 
theory for RR states at v = k/(kM + 2) can be described 
as RRfe m = U(l)fe ( jvf x Pfj. where U(i)k,M is due to the 
electric charge and Pfj. represents the Zj.-parafermion the- 
ory JHEl. The U(l) k,M part °f this theory is a simple 
Abelian contribution, essentially labeled by integral multi- 
ples of the charge/flux unit ( km+i ' If)' wnere — e is the 
electron charge and — = <I>o is the magnetic flux quan- 
tum (in units H = c = 1). The fusion rules for these 
labels are just addition of charge/flux and the S-matrix is 
S na n b = e'"-" l W5/ v / p+2). 

The Pffe part of this theory requires more explanation (for 
a discussion of its braiding, see 13 811 '). Essentially, we use 
that the theory is equivalent to the coset SU(2) fc /U(l) fc . As 
a consequence, the Z^-parafermion sector's anyonic charge 
can be labeled by the corresponding CFT fields where 
A G {0, 1, . . . , fc} and A 6 Z, subject to the identifications 



ff> A 

V A+2fc 



^A-fc and the restriction A + A = (mod2) 
(giving a total of hk (k + 1) fields). The fusion rules for this 
sector are (as for the CFT fields): 



a{Ao+Ai,,2fc-A £ ,-A i ,} 



$ A„ x $ A 6 = 



E 



4> 



A 

A +A 5 



A=|A a -A b | 

Their quantum dimensions are 

'(A + 1)tt s 



d$A = sin 



fc + 2 



fc + 2 



Special fields in this theory are the vacuum 1 

parafermions ipi = the primary fields 07 

Zfc-neutral excitations ej = <fe^ , where 1=1, 

3 = 1, • • 
for Pffc is 



(8) 



(9) 

*g, the Z fe - 
$5, and the 



-/' 
,k 



1 and 



\_(k — 1) /2J. From ijfji we find that the S-matrix 
sin ^ 



(Ag+l)(Aj,+l)7 

Ai+2 



„— iA a \ b n/k 



(10) 



The U(l)fc,M and Pffc sectors combine so that the any- 
onic charges in the RRfc j\/ theory are (defining a shorthand) 

^(fclTI-t 1 '^")' where we have ™ eZ and A « G 
{0, 1, ... , k} such that A„ + n = (mod2). Quasiholes carry 

charge ii = ( h .^ e ,^ 7 , ff , 01 The S-matrix for RR anyons 



is obtained by multiplying the S-matrix elements of the two 
sectors for anyons a and b, and renormalizing by some overall 
constant c (which we will not need explicitly): 



S ab = csin 



/ (A na + l)(A» b + l)7r 
V fc + 2 



e -in a n b hM + 2 , (11) 



Since the tunneling current is dominated by quasiholes |4 
we only need the monodromy matrix elements 



M 



(A„ + 1)7T 

fc+2 



li,A„ 



(12) 



We note that RR24 is the MR state, and we can easily check 
that this exactly matches the results of 

We now turn to the RR^i theory for v — 12/5. The Pf3 
theory has six fields: 1, ipi, ip2, which have quantum di- 
mension 1, and cti, 02, e, which have quantum dimension 
(the golden ratio). The total quantum 



= 2 cos ■ 



1W5 



dimension is D = y^3 (4> + 2) and the S-matrix is 



S =D 



1 

e 3 



1 

e 3 
e 3 



e 3 
-1 



<pe 3 

e~ l "3 
-1 



-1 
-1 
-1 



where the columns and rows are in the order: 1, ipi, tp2> ci, 
<72, e. Quasiholes in the RR3,i theory have anyonic charge 
(f j |f i o'l)- It is useful to consider a Bratteli diagram (which 
has periodicity 6 in n) to keep track of the allowed Pf3 charge 
for a corresponding value of n: 
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The longitudinal conductivity in the interferometry experi- 
ment will be 



<Jxx oc |ii| 2 +|£2| 2 +2 \tit 2 



l(-r 2 )%os(/3 + nf) 



(13) 



where = 1 if the n quasihole composite on the antidot has 
Pf3 charge with quantum dimension <p (i.e. <7i, 02, or e) and 
= if the composite has quantum dimension 1 (i.e. 1, ipi, 
or if>2)- Thus, depending on the total Pf3 charge on the antidot, 
one of two possible conductivity values will be observed. The 
Pf 3 charge may then be determined by varying j3 (using the 
side gate) to measure the interference fringe amplitude, which 
is suppressed by a factor of $r 2 w .38 when = 1, This 
behavior indicates the presence of non- Abelian statistics, and 
distinguishes this state from other proposals for the same fill- 
ing fraction (e.g. composite fermions). To describe v = 12/5, 
we apply a particle-hole transformation to RRs^, replacing S 
with S (hence, M with M), which changes the sign in front 
of n in Eq. Jl 3i 1 46] . In more general scenarios where com- 
posites of quasiholes/quasielectrons may be used instead of 
the single tunneling quasiholes, identical behavior (up to the 
phase) will be exhibited by quantum dimension <f> composites, 
while quantum dimension 1 composites exhibit a single un- 
suppressed interference pattern at any n. 

We conclude with a few remarks explaining that despite 
the relatively simple nature of these interferometry experi- 
ments, they provide a surprisingly large amount of informa- 
tion about the system being probed. This is because the 



4 



experiments essentially measure the S-matrix of the TQFT 
that describes the system. The S-matrix fully determines 
the fusion rules through the Verlinde formula 1 3 



N c , = 

ab 

J2xec SaxSbxScx I Sf x . Additionally, a theorem known as 
"Ocneanu rigidity" states that, given a set of fusion rules, there 
are only finitely many corresponding TQFTs with these rules 
14011 . In other words, knowledge of the S-matrix is sufficient 
to pin down the topological order of the state to a finite num- 
ber of possibilities. Clearly, it may be difficult to measure all 
elements of the S-matrix by the methods described here. It 
appears particularly challenging to invoke tunneling of any- 
onic charges different from that of the quasihole, though one 
may speculate on techniques that may eventually prove suc- 
cessful, such as resonant effects with intermediate antidots of 
tunable geometry and capacity on the tunneling arms. Still, 
the S-matrix has many special properties and so even a partial 
measurement of fairly low accuracy may be sufficient to deter- 
mine it. In addition to Eq. 0, any S-matrix must satisfy a set 
of constraints coming from the Verlinde formula and the fact 
that the fusion coefficients are integers. Also, the first row of 
the S-matrix must be real and positive, because of its relation 
to the quantum dimensions (and in fact, all elements are num- 
bers with special algebraic properties). Finally, given an S- 
matrix, there must be a diagonal matrix T which together with 
S generates a representation of the modular group SL(2, Z), 
implying that (ST) 3 = S 2 . For any fixed number of parti- 
cle types, only finitely many different S-matrices are known 
(and it is conjectured that only finitely many exist). Hence, 
once some the S-matrix elements and the number of different 
charges are known from measurements, one may look at the 
finite list of known S-matrices and hope to identify one that 
matches. In conclusion, for any two-dimensional system, in- 
terference experiments as described here can in principle de- 
termine the fusion rules and even a finite set of TQFTs, one of 
which will fully describe the topological order. 
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[43] In the weak tunneling regime, the tunneling current I oc y 2s_1 
where s is the scaling dimension/topological spin of the corre- 
sponding fields 1 41]. It follows that the dominant contribution 
is from the field with lowest scaling dimension, which in FQH 
systems is the quasihole. 

[44] We expect the states of the antidot with the same electric charge 
but different topological spins to be non-degenerate (with en- 
ergy difference scaling as L^ 1 for an antidot of size L) and to 
have different charge distributions, ruling out the possibility of 



their superposition over extended time. 
[45] As previously noted, the particle with lowest topological spin 
dominates tunneling. For RRfc.M particles with charge A„, the 
topological spin is given by s An = 2k( f M+2) + - 

where the field identifications must be used to write — $^ 
such that —A < A < A in order to apply this formula. One 
may check that all other particles have greater spin than the 
quasihole value, a u = 2^+2)^X2) • 
[46] We thank S. B. Chung for bringing this sign change to our at- 
tention. 



